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SYSTEM AND METHOD OF POINT MATCHING MEASURED 
POSITIONS TO TEMPLATE POSITIONS 

CROSS-REFERENCE TO RELATED APPLICATIONS 

The present application claims the benefit of United States provisional application 
5 Serial No. 60/454,541, filed March 14, 2003, entitled "A METHOD OF FAST, ROBUST, 

RIGID POINT MATCHING OF PROBES TO TEMPLATE POSITIONS," the disclosure of 
which is hereby incorporated herein by reference in its entirety. 

FIELD OF THE INVENTION 

Aspects of the present invention relate generally to the field of matching sets of 
10 points between a measured set and a template set, and more particularly to a system and 

method providing fast, robust, rigid point matching of measured positions to template 
positions. 

BACKGROUND OF THE INVENTION 

Matching sets of points between a measured set and a template set is a common 
15 problem in astronomy, biomedical imaging, animation, precision instrument testing and 

calibration, and other technologies. One such technology is probe card analysis, in which 
planarity and alignment of a plurality of probes arranged in an array may be tested or 
calibrated through comparison with a known location of one or more fiducial marks, for 
example, disposed on or integrated into the structure of a fiducial plate. In particular, some 
.20 probe card analysis systems employ optical or image data acquired by a camera or other 

imaging apparatus; in some cases, the imaging apparatus may be disposed on the opposite 
side of a substantially transparent fiducial plate from the probe array associated with the 
probe card. These acquired data are representative of both probe tips and fiducial marks, and 
data processing may be employed which seek to identify position and orientation of probe 
25 tips relative to each other through comparison with the known position and orientation of 

one or more fiducial marks. Alternatively, the position and orientation of imaged probe tips 
may be compared to ideal position and orientation data maintained in the template or 
reference data set. 

In the case of a probe card analysis system, image data representing the probe tips 
30 may generally be considered "measured" data or points; in some applications, the "template" 

data may represent ideal or reference points (derived, for example, from the design 
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specifications of the probe card or an industry standard) independent of imaged fiducial 
marks. Similarly, many of the foregoing technologies, among others, employ comparative 
metrology relying upon measured data sets and template, or reference, data sets; in many 
instances, point matching techniques are employed in an attempt to match data representing 
5 non-ideal measured points (which include errors, missing points, and erroneous points) to 

data representing known ideal, template, or reference points. In practice, one example of an 
application in which point matching techniques may have utility is in comparing an image 
acquired by a telescope to known locations of celestial objects as recorded on a star chart. 
While some data processing components attempt point matching operations to fit the 
10 measured points to a known reference frame, currently implemented systems and methods 

are deficient in that traditional methodologies do not accommodate fast, robust, rigid point 
matching of measured positions to template positions in an integrated and elegant solution. 

SUMMARY 

Aspects of the present invention overcome the foregoing and other shortcomings of 

15 conventional technology, providing a system and method of point matching measured 

positions to template, or reference, positions for various applications. Exemplary systems 
and methods may generally comprise creating force field vectors and moments which are 
operative to perturb measured point locations into alignment with template point locations. 
Progressive reduction of the force radius may produce results that match very well for points 

20 that are highly correlated with the template pattern. 

In accordance with some embodiments, for example, a method of point matching 
measured points to template points comprises: acquiring measured data representing a set of 
measured point locations; comparing the set of measured point locations to template data 
representing a set of template point locations; defining force field vectors operative to 

25 perturb the measured point locations into alignment with the template point locations; and 

responsive to the defining, matching measured point locations to template point locations. 

Depending upon the application or particular system requirements, the matching 
may comprise utilizing a many-on-many point matching algorithm or a one-to-one point 
matching algorithm. A method employing one such algorithm may additionally comprise 

30 determining offsets, rotations, and position errors between the measured point locations and 

the template point locations. 
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In some implementations, the defining may comprise creating the force field vectors 
to act over a prescribed range. In one the creating may comprise, for distances greater than 
the prescribed range, dissipating a magnitude of the force field vectors with increasing 
distance; additionally, the creating may comprise, for distances shorter than the prescribed 
5 range, increasing a magnitude of the force field vectors with increasing distance. 

Some methods of point matching measured points to template points may further 
comprise selectively repeating the comparing, the defining, and the matching operations set 
forth above; additionally, the selectively repeating may comprise selectively decrementing 
the prescribed range at each successive iteration of the defining. 

10 In accordance with another aspect of the subject matter described herein, a computer 

readable medium encoded with data and instructions for point matching measured points to 
expected points may cause an apparatus executing the instructions to: acquire measured data 
representing a set of measured point locations; compare the set of measured point locations 
to reference data representing a set of expected point locations; define force field vectors and 

15 moment arms operative to perturb the measured point locations into alignment with the 

expected point locations; and selectively repeat, until predetermined convergence criteria 
have been satisfied: comparing, to the reference data, the measured point locations perturbed 
by the force field vectors and the moment arms; and redefining the force field vectors and 
the moment arms responsive to the comparing. 

20 In some applications, such a computer readable medium may further cause an 

apparatus executing the instructions to compare the set of measured point locations to the set 
of expected point locations utilizing a many-on-many or a one-to-one point matching 
algorithm. Such instructions may further cause an apparatus to determine offsets, rotations, 
and position errors between the measured point locations and the expected point locations. 

25 As set forth above, the force field vectors may be created to act over a prescribed 

range; in that regard, a computer readable medium may further cause an apparatus executing 
the instructions to decrement the prescribed range at each successive iteration of the 
redefining. 

In accordance with the principles set forth in detail below, an exemplary method of 
30 measuring probe locations in a probe card analyzer system may comprise: acquiring 

measured data representing a set of probe point locations; comparing the set of probe point 
locations to reference data representing a set of expected point locations; responsive to the 
comparing, defining force field vectors and moment arms operative to perturb the probe 
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point locations into alignment with the expected point locations; and selectively repeating: 
computing modified probe point locations representative of the probe point locations 
perturbed by the force field vectors and the moment arms; and redefining the force field 
vectors and the moment arms responsive to the computing; until predetermined convergence 
5 criteria have been satisfied. The acquiring may comprise utilizing an imaging apparatus. 

The computing may comprise determining offsets, rotations, and position errors 
between the probe point locations and the expected point locations. A method of measuring 
probe locations in a probe card analyzer system is contemplated wherein the defining and the 
redefining comprise creating the force field vectors to act over a prescribed range; in some 
10 applications, the redefining comprises decrementing the prescribed range at each successive 

iteration of the redefining. 

The foregoing and other aspects of the disclosed embodiments will be more fully 
understood through examination of the following detailed description thereof in conjunction 
with the drawing figures. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a simplified diagram depicting one embodiment of a general case force 
field definition and illustrating active forces for two measured points. 

FIG. 2 is a simplified diagram illustrating one embodiment of a special case force 
field involving a one-to-one point matching scenario. 
20 FIGS. 3-8 illustrate results achieved by one embodiment of a point matching method 

and respectively depict measured points becoming progressively more aligned with template 
points at the conclusion of successive major iteration loops. 

FIGS. 9-14 illustrate results achieved by one embodiment of a point matching 
method and respectively depict measured points becoming progressively more aligned with 
25 template points at the conclusion of successive major iteration loops. 

DETAILED DESCRIPTION 

By way of background, it is noted that the exemplary algorithm set forth below may 
provide or enable a fast, robust method for performing rigid point matching between 
measured and template data sets, and for determining the rigid transformation required to 
30 bring the two sets into best alignment. In that regard, an exemplary point matching strategy 

matches measured locations of points to expected locations of points. The exemplary 
algorithm is robust to outliers and missing or extra data in the measured data set, and may 
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generally perform a many-on-many point matching operation where the correspondence 
between points is unknown. In accordance with some simplified embodiments, an 
exemplary algorithm may also be modified to perform one-to-one point matching when a 
one-to-one point correspondence is known. Those of skill in the art will appreciate that 
5 other algorithmic approaches or recursive computational procedures may provide similar 

functionality and may readily be implemented in accordance with the principles described 
herein. The present disclosure is not intended to be limited to any particular equation or set 
of equations having utility in achieving the desired results. 

The point matching functionality set forth herein may employ both measured data 

10 sets and template, or reference, data sets. Specifically, a system and method configured and 

operative in accordance with the present disclosure match data representing possibly non- 
ideal, measured points (i.e., "measured" data may include errors, missing points, and 
erroneous points) to data representing known ideal points. In that regard, the terms 
"template," "ideal," and "reference" as used herein are generally intended to refer to a set of 

15 known or given data points against which "measured" or acquired data points are measured; 

in this context, it will be appreciated that the term "template" is not necessarily intended to 
connote physical structure (such as a fiducial plate or a printed star chart, for example), 
although it may in some applications. 

As set forth in more detail below, an heuristic approach operative in accordance with 

20 the present disclosure may create an artificial, non-linear, attractive "force field" between all 

possible point pairs in the general case. FIG. 1 is a simplified diagram depicting one 
embodiment of a general case force field definition and illustrating active forces for two 
measured points. As represented in FIG. 1, the "O" symbols denote template points which 
are fixed, and the "X" symbols denote measured points sought to be brought into optimal 

25 alignment with the template points. In accordance with one aspect of the present invention, 

each measured point "X" may be considered to be drawn, or "perturbed," toward all 
template points by the force field. In addition to overall translation and rotation, errors in the 
placement of the measured points may also exist. Such errors may be introduced, for 
example, by resolution characteristics, acquisition techniques, or other limitations associated 

30 with the imaging apparatus or other hardware employed to acquire the measured data set. 

Additionally, FIG. 1 illustrates a condition in which a point is missing from the measured 
data set, le. 9 in the exemplary FIG. 1 embodiment, only four measured points are present, 
whereas five points are present in the template set. 
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In practice, the one-to-one matching situation may represent a special case; such a 
special one-to-one point matching case may utilize what are referred to herein as "attractive 
forces" only between corresponding points, as generally illustrated in FIG. 2. In particular, 
FIG. 2 is a simplified diagram illustrating one embodiment of a special case force field 
5 involving a one-to-one point matching scenario. In the FIG. 2 embodiment, a known one-to- 

one correspondence exists between template points and measured points, Le. 9 each respective 
measured point is known to correspond to a respective template point. 

When summed, resultant attractive forces and moments (embodiments of which are 
defined by way of example below) may be used to bring the measured point set into 

10 alignment with the template point set. The force field vectors (e.g., represented by the 

dashed arrows in FIG. 2) may be created to act over a prescribed range, r. For distances 
greater than r, the attractive forces may be designed to dissipate with increasing distance; for 
distances less than r, however, the attractive forces may be designed to increase with 
increasing distance. In accordance with some embodiments, the foregoing force field 

15 properties may allow an algorithm such as set forth herein, or a similar data processing 

technique, robustly to handle, model, or otherwise to accommodate outlying, missing, or 
extra data points, while tightly matching points with high correlation between data sets. 

In some instances, the solution process may be configured and operative iteratively 
to apply rigid shift and rotation to the measured point set in the direction of the overall forces 

20 and moments. In that regard, an exemplary solution process may initially create or model a 

force field having a large range of influence, and may additionally execute iterations until 
the measured and template point centroids and principle axes align. The range of influence 
of such a modeled force field may then be made progressively smaller so as to fit points with 
high correlation between sets even more tightly. 

25 In accordance with some exemplary embodiments, the algorithm described herein 

may be a general algorithm for points in n-space. Without loss of generality, however, and 
for the sake of clarity, the following description generally makes reference to the two- 
dimensional case, z.e., wherein n = 2. 

Many-on-Many Point Matching . 

30 It will be appreciated that a given template set may comprise a set of points having 

coordinates (x^ , y ti ), where (i = 1, 2, . . . , n pt ). Similarly, a measured point set may generally 
comprise points having coordinates (x^ , y m j), where ( j = 1, 2, . . . , n pm ). For the set of 
measured points, an exemplary algorithm may calculate a centroid (x™ , y m ), a mean central 
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radius, r m , and a radius of gyration, p m , for example, in accordance with the following 
relationships: 
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fw = 



Tm = 



U P» ■/-!■ 

I n Pm 

^ I >(y«. 



,(2) 
(3) 



Pw = 



In one embodiment, the variables w bb and h bb may define a width and height, 
respectively, of a bounding box of the set of measured points. The maximum range of 
influence of the force field, also referred to as the maximum force radius, in this case may 
then be defined as 



(5) 



where r fi i mit is a user prescribed (or otherwise preselected or predetermined) limit on the 
maximum allowable range of influence. The range of influence may then be progressively 
reined in or incrementally diminished (i.e., selectively decremented) over major iteration 
steps from an initial value, r^, to a final value, r^. If the variable k, defines the reduction 
factor in the force field range of influence from one major iteration loop to the next, then the 
number of force radius reductions required to reach r^ may be expressed as: 



M, P int , ^ + 1 



(6) 



Exemplary Force Field Definition 
Where the current major iteration step index is denoted by a variable, m, the force 
radius at this current major iteration step may then be set in accordance with Equation (7) set 
forth below. 
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Between every point, i, in the template set and every point, j, in the measured set, 
normalized distance vector components between these points may be expressed as follows: 

(8) 

a^ ^zM , ; (9) 

It is noted that the foregoing distance vector components are normalized by the force 
radius r f . This normalization has the effect of allowing distances shorter than the force 
radius to produce normalized distances of less than one, and allowing distances longer than 
the force radius to produce normalized distances greater than one. As set forth briefly above, 
such normalization may provide the force field with properties and characteristics useful for 
handling outlying and missing points from the measured point data set. The normalized 
distance between points i and j may then be described as indicated in Equation (10). 

^ 4 v^ ; (1D> 

It is also convenient to define moment arms for each point-to-point force vector. In 
an exemplary implementation, each force vector represents the force applied to a particular 
measured data point and has a moment arm which may be computed about the centroid of 
the measured point array (Le. 9 all the points in the measured data set). In accordance with 
this embodiment, 

(id 

The normalized point-to-point force vector components of the force on measured 
point j exerted by template point i may then be defined as follows: 



f«i= , \ 2 (13) 



(14), 



Wm : 

r"ii=ixifn-h^> (15) 

For normalized distances that are large compared to one, the force vector components 
diminish quickly as 1/r 3 , allowing robust rejection of missing data points or outliers. For 



8 



normalized distances that are small compared to one, the force vector components become 
linearly increasing with increasing distance, enabling smooth convergence properties as the 
pattern match process progresses to a final solution. 

In some embodiments, the foregoing point-to-point forces may be aggregated into an 
overall average force and moment applied to the measured point cluster. 

it tiff ttpm 

- : (16) 

\ n pt n vm - ..- t, ' 

. n pi n Pm i=i/=i :> . . • ■ ^l,^;:; 

«=— X Tm (is) 

Iterative Solution Process 
At each minor iteration step, aggregate forces F x , F y , and M may be calculated, for 
example, as set forth above or in a similar manner. In this context, the essence of a "minor" 
iteration step is to perturb the measured point cluster in the direction of the applied forces 
with the aim of eventually bringing the measured points making up the cluster into optimal 
alignment with the template points. Translation and rotation perturbations about the centroid 
of the measured point array may be expressed as follows 

|a = %$T : :,, (21) 

where k x , k y , and ke represent the convergence gains. These gains may be selected to 
provide rapid, stable convergence, and may be application or implementation specific. In 
one embodiment, optimal gain selection may be as expressed in Equations (22) through (24). 

'' ; fk x = : tf {72) 

^ = 7 .-.p.- '' . 11!|(23> : : 

WM ^ '11 (24) 



This minor iteration loop may be selectively repeated, or may automatically 
continue, for example, either until a convergence criterion (or more than one convergence 
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criteria) has been satisfied, or until a predetermined or preselected number of iterations have 
been performed. One exemplary convergence criterion may be defined as the final 
prescribed force radius: 

= ^ r v ' : ; v ' : . ^ (25) 

In some embodiments, the test for convergence may be satisfied if all of the 
following conditions are met. 

6 x <t (26) 
■ &y<t (27) 
V™ < t (28) 

Once convergence has been achieved, or following a maximum number of iterations, 
the minor iteration loop set forth above may terminate, and the next major iteration loop may 
begin. In this context, a "major" iteration loop may function progressively to reduce the 
force radius as set forth above. As the force radius decreases, the solution becomes more 
attuned to matching measured points that align most closely with the template. At the 
completion of the defined major iteration loops, i.e., a predetermined or preselected number 
of major iteration loop operations, the measured point set may be optimally aligned with the 
template point set. 

FIGS. 3-8 illustrate results achieved by one embodiment of a point matching method 
and respectively depict measured points becoming progressively more aligned with template 
points at the conclusion of successive major iteration loops. Similarly, FIGS. 9-14 illustrate 
results achieved by one embodiment of a point matching method and respectively depict 
measured points becoming progressively more aligned with template points at the conclusion 
of successive major iteration loops. 

Results of the iterative process depicted in FIGS. 3-8 are representative of a simple 
case involving five template points, and four measured points (such as the situation depicted 
in FIG. 1, for example). Results of the iterative process depicted in FIGS. 9-14, on the other 
hand, are representative of a more complicated point matching example. As the progressive 
results indicate in FIGS. 8 and 14, the measured points may include errors which prevent 
them from perfectly aligning with the template points; some possible sources of such errors 
are addressed above with specific reference to FIG. 1. Note that in both examples, however, 
the effect of progressive reduction of the force radius may produce results that match very 
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well for points that are highly correlated with the template pattern-there is little evidence of 
trading large errors on low correlation points with small errors on high correlation points. 



In accordance with one variation of the general algorithm set forth above, 
knowledge of correspondence between measured data points and template data points may 
be assumed. In such a case, the attractive force field may be set up only between 
corresponding points in the template data set and the measured data set. This simplification 
greatly reduces the computational complexity of the algorithmic approach (0(n 2 ) — > O(n)), 
and consequently allows a solution to converge more quickly than in the general case where 
correspondence is unknown or assumed not to exist. 

The preliminary calculations described in Equations (1) through (6) may remain 
unchanged for a one-to-one variation of the algorithm. The force definition parameters in 
Equations (7) through (15) may be modified for one-to-one point matching by setting i = j. 
Then the total average force between the two data sets may be expressed as follows: 



The iterative solution process may proceed substantially as set forth above with 
reference to the many-on-many point matching scheme. 

Those of skill in the art will appreciate that the presence of unique, non-repeating 
features may affect identification of a global optimum solution. If such unique features are 
not present, a point matching algorithm such as described herein may converge to a local 
minimum. The stronger the presence of unique features, the more likely the algorithm is to 
identify a global optimum. 

Aspects of the present invention have been illustrated and described in detail with 
reference to particular embodiments by way of example only, and not by way of limitation. 
It will be appreciated that various modifications and alterations may be made to the 



One-to-One Point Matching 
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exemplary embodiments without departing from the scope and contemplation of the present 
disclosure. It is intended, therefore, that the invention be considered as limited only by the 
scope of the appended claims 
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